Elliptic flow of deuterons in relativistic heavy-ion collisions 
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Abstract 

Using a dynamical model based on the A^A'^ — > dn, NNN dN, and NNir — > dn reactions 
and measured proton and pion transverse momentum spectra and elliptic flows, we study the 
. production of deuterons and their elliptic flow in heavy ion collisions at RHIC. The results are 

I compared with those from the coalescence model. The deviation of deuteron elliptic flow from 

Q ' the constituent nucleon number scaling expected from the coalescence model and the comparison 

Q ■ with the experimental data are discussed in connection to the allowed nucleon phase space in these 

\^ • reactions. 
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I. INTRODUCTION 



The elliptic flow of particles in heavy ion collisions is a measure of the strength of the 
second Fourier coefficient in the azimuthal angle distribution of particle transverse momen- 
tum relative to the reaction plane [1, 2]. Significant information on the properties of the 
hot dense matter formed during the initial stage of heavy ion collisions has been obtained 
from the study of elliptic flow [3-5]. In heavy ion collisions at the Relativistic Heavy Ion 
Collider (RHIC) [6, 7], measured elliptic flow at low transverse momentum shows a mass 
ordering, i.e., the strength of elliptic flow of identifled hadrons decreases with increasing 
hadron masses, and this has been well described by ideal hydrodynamics [8] as well as by 
the transport model [9]. Furthermore, the elliptic flows of identified hadrons, particularly 
at intermediate transverse momenta, are seen to follow a constituent quark number scal- 
ing, i.e., the dependence of hadron elliptic flows on hadron transverse momentum becomes 
similar if both are divided by the number of constituent quarks in a hadron. This scaling 
behavior of hadron elliptic flows is consistent with the predictions of the quark coalescence 
model for hadron production from produced quark-gluon plasma in relativistic heavy ion 
collisions [10-14]. 

Recently, the elliptic flow of deuteron has also been measured at RHIC [15-17]. These 
measurements show that the deuteron elliptic flow seems to scale with its constituent nu- 
cleon number, implying that the quark number scaling of elliptic flows holds not only for 
hadrons but also for the deuteron. Since deuteron production in heavy-ion collisions could 
be described by the coalescence of protons and neutrons at freeze out [18], the observed 
nucleon number scaling of the deuteron elliptic flow is thus not surprising. In the simplest 
coalescence model [13, 14] which involves only comoving particles, the deuteron yield in 
momentum space is proportional to the product of the proton and neutron densities at half 
the momentum of produced deuteron, i.e., Pp = Pn = Pd/'2, and the deuteron elliptic flow 
would satisfy exactly the nucleon number scaling and thus the quark number scaling as well. 
This strong momentum constraint is relaxed in the more general coalescence model [18] that 
takes into account the nucleon momentum spread in the deuteron. As in the case of hadron 
production from quark coalescence [10], the more general coalescence model would lead to 
a small deviation of the scaled deuteron elliptic flow from that of nucleons. 

Although the momentum conservation is maintained in the coalescence model, the energy 
conservation is not satisfled. Some doubt has thus been raised on the applicability of the 
coalescence model, especially in low transverse momentum region [19, 20]. Investigating 
the consequences arising from the energy conservation condition is therefore required to 
test the coalescence model and to understand the constituent number scaling of elliptic 
flows in low transverse momentum region. In this work, this question will be addressed by 
studying deuteron production and elliptic flow in heavy-ion collisions based on a dynamical 
approach that satisfles both energy and momentum conservations. Speciflcally, deuteron 
production in the present study will be treated through NN dn, NNN — > dN, and 
NNtt — > dn reactions. Since this approach is based on physical scattering processes for 
deuteron production, the energy-momentum conservation is always satisfied. By comparing 
the physically allowed nucleon phase space for deuteron production with that involved in the 
coalescence model, we can study the consequences of the energy conservation condition. In a 
realistic study, these reactions need to be included in the hadronic stage of relativistic heavy 
ion collisions either via a microscopic transport model [21] or a schematic kinetic model [22] 
to take into account all possible deuteron production and dissociation processes. Since the 
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elliptic flow of particles during the hadronic evolution does not change much as the initial 
spatial asymmetry through which the elliptic flow is generated has decreased signiflcantly 
during the evolution of the partonic stage [9], we can therefore calculate the deuteron yield 
and elliptic flow based on the deuteron production rate from nucleons at freeze out. 

This paper is organized as follows. In Section II, we flrst give a brief description of the 
coalescence model that will be used to compare with the results from the dynamical model 
employed in the present study for deuteron production. We then consider in Section III 
and Section IV deuteron production from the two-body reaction NN — > dn and three-body 
reactions NNN — > dN and A'"A'"7r — > dn, respectively, with their cross sections obtained from 
a hadronic model that is based on empirically determined parameters [23]. In Section V, 
nucleon and pion transverse momentum spectra and elliptic flows, which are needed in 
both the coalescence model and the dynamical model for studying deuteron production, are 
discussed. Results on deuteron transverse momentum spectrum and elliptic flow are then 
compared with available experimental data from RHIC as well as those from the coalescence 
model in Section VI, where we further compare the nucleon phase space involved in our 
dynamical model with that in the coalescence model. Section VII contains the summary 
and discussions. Details on the reaction amplitudes for both two-body and three-body 
reactions are given in Appendixes A and B, respectively. 



II. DEUTERON PRODUCTION IN THE COALESCENCE MODEL 

In the coalescence model, the probability for the production of a bound composite particle 
from a many-particle system is determined by the overlap of the wave functions of coalescing 
particles with the internal wave function of the composite particle [24]. Assuming that 
nucleons are uniformly distributed in space, the momentum spectrum of deuterons formed 
from the coalescence of nucleons is given by 
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In the above, the factor 3/4 is the probability for a proton and a neutron to form a spin triplet 
state like that in a deuteron; fp{Pi) and fn{P2) ^"^^-i respectively, the proton and neutron 
momentum distributions in the hadronic matter of volume V at freeze out, and they are 
related to the proton number Np and neutron number Nn by V J -^;^fp,n{p) = Np^n', is 
the deuteron wave function in the momentum space with p[ and denoting, respectively, 
the momenta of the proton and neutron in the deuteron rest frame. In the more general case 
of non-uniform nucleon distribution and/or with collective flow, additional spatial integrals 
would appear in Eq. (1), and the deuteron momentum space wave function is replaced by 
its Wigner function which describes both the relative momentum and spatial distributions 
of the two nucleons in a deuteron [25, 26]. 

For the deuteron wave function, we use the one given by Hulthen. In the momentum 
space, it is given by 

with = 0.23 fm~^ and f3^ = 1.61 fm~^ [26] and is normalized as / d'^k\'^d{k)\'^ = 1. It 
should be mentioned that this more general coalescence model still does not respect the 
energy conservation condition. 
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In the naive coalescence model, the deuteron wave function in the momentum space is 
taken to be a delta function. In this case, the transverse momenta prp^ and ^^"2 of ^^e 
nucleons forming the deuteron are restricted not only in magnitude but also in direction 
so that prpi = Pt2 — PtJ'^ with pj, ^ being the deuteron transverse momentum [13, 14]. 
As a result, the scaling of the nucleon elliptic flow U2,jv and the deuteron elliptic flow V2 ^ 
according to their constituent nucleon numbers is exact, i.e., V2^J(Pt) — '^'^2,NiPT/'^) [1^, 2?]. 

III. DEUTERON PRODUCTION FROM TWO-BODY REACTIONS 

For deuteron production from two-body nucleon-nucleon reactions, the dominant reac- 
tion is NN — > dn due to the small pion mass and the strong pion-nucleon coupling. Indeed, 
experimental observations have shown that the cross section for the NN dn reaction 
is much larger than the cross sections for deuteron production from other two-body reac- 
tions [23]. In contrast to the coalescence model, in which the deuteron is formed from a 
proton and a neutron in the hadronic matter, not only the pn reaction but also the pp and 
nn reactions can produce a deuteron in the NN — > dir reaction as the pion, which has 
isospin one, can carry away the mismatched isospin. 

The production rate of deuterons with momentum p^ from the NN — > dn reaction can 
be written as 

dTZd _ 1 /■ -A d^Pi d^p^ f ( -)] 

X \M{NN ^ d7r)r(27r)^5(^) {flP^ ' ' Pd^^ ' (3) 

where fN{p) and fnip-n) are nucleon and pion momentum distribution functions with the 
latter normalized to the pion number in the hadronic matter. We have assumed here that 
the proton and neutron densities in the hadronic matter are the same for simplicity. The 
average of initial spin and isospin states and sum over final spin and isospin states are 
included in the squared transition amplitude \M.{NN (i7r)p. Details on the amplitude 
for the NN — > (ivr reaction are given in Appendix A. Compared with the coalescence model, 
Eq. (3) shows that both energy and momentum are always conserved due to the produced 
pion. This allows nucleons with higher momentum to produce a low-momentum deuteron 
since the pion can carry off the excess energy of the two nucleons. As a result, there is no a 
priori restriction on the momenta such as pj,^^ ~ 2py as in the coalescence model. Instead, 
the nucleon distribution and the energy-momentum conservation condition determine the 
most plausible configurations of nucleon momentum for producing the deuteron. 



IV. DEUTERON PRODUCTION FROM THREE-BODY REACTIONS 

Deuterons can also be produced from three-body reactions such as NNN dN and 
NNir — > dir. The reaction NNN — > dN is known to dominate deuteron production in low- 
energy heavy ion collisions [28, 29]. Because of the large abundance of pions in the produced 
hadronic matter in heavy ion collisions at RHIC, the reaction NNn — > dn becomes also 
possible. 
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A. NNN dN reaction 



For the NNN — > dN reaction, the production rate for deuterons of momentum is 
given by 

dTZd _ 1 f -|\ d^Pi , , ^ d^PN n / / m 
" {2nr2Ea J (27r)32£;/^^^^^ (27r)32£;^, ~ ^^^^^^^ 

X |A1(iVArAr ^ dAr)|2(27r)^5(^) (j2p,-p^-p}j . (4) 

The presence of the third nucleon in the initial state leads to a very different nucleon 
phase space in deuteron formation, as the third (spectator) nucleon can ensure the energy- 
momentum conservation. Since the nucleon mass is much larger than the pion mass, the 
effect of the spectator nucleon is different from that of the pion produced in the A^A^ — > 
dn reaction and thus could affect the elliptic flow of produced deuteron differently. The 
amplitude for the NNN — > dN reaction can be obtained from its inverse reaction dN — > 
NNN, which has been studied in Refs. [29-31] based on the full NN elastic scattering 
amplitude. Instead of such a complete study, wc approximate the amplitude for the reaction 
NNN — > dN by extending the amplitude wc derived in the previous section for the NN — >■ 
dTT reaction to an off-shell pion that is absorbed by the third spectator nucleon. Details on 
the modehng of the amplitude A4{NNN — > dN) are given in Appendix B. 



B. NNtt d-K reaction 



The same method introduced in the above for the NNN — > dN reaction can be used 
for the amplitude of the NNn dn reaction, i.e., the subprocess A^A^ dn is replaced 
by A^A^ — > dp with the off-shell rho meson absorbed by the spectator pion. Details on 
the modehng of the amplitude M.{NNt^ dn) are also given in Appendix B. In terms 
of the nucleon and pion momentum distributions the production rate for deuterons with 
momentum from this reaction is given by 

d'R'd _ 1 /• d^Pi d^p^ £'p^ 

d^Pd (27r)32Ed J A| (27r)32E/'^^^^^ {2'n)nE/^^'^^' (27r)32E,^ ^ ^""^^^'^ 

X \M{NNn d7r)\\27r)'S^'^ (j^p, -p^-p^ . (5) 



V. NUCLEON AND PION TRANSVERSE MOMENTUM DISTRIBUTIONS 



Both the coalescence model and the dynamical model use as input the nucleon momentum 
distribution at freeze out. Moreover, the pion momentum distribution at freeze out is also 
needed in the dynamical model. For heavy ion collisions at RHIC, both nucleon and pion 
rapidity distributions are approximately uniform in the mid-rapidity region, i.e., \y\ < 0.5, 
which we will assume in the present study. In this case, only their transverse momentum 
distributions are needed. In the following, we discuss how these distributions are obtained. 
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FIG. 1: (Color online) (a) prp spectrum of protons. The solid line is half of the parameterized 
nucleon distribution [Eq.(6)]. Experimental data are from Ref. [32] for the minimum-bias pro- 
ton spectrum (filled square), (b) Nucleon differential elliptic flow V2{pj')- The solid line is the 
parameterized one [Eq. (7)] and filled squares are experimental data from Ref. [16]. 



A. Nucleon spectrum and elliptic flow 

For the nucleon pj, distribution fNiPr)^ "we take it to have the form 

fNiPr) = 7ivexp (-m^/Teff) [1 + v^^NiPr) cos(20)], (6) 

where rrirp = a/ Af ^ + with the nucleon mass and being the azimuthal angle of 
the nucleon transverse momentum with respect to the reaction plane. We have introduced 
a dimensionless parameter 7^ to compensate for the effects caused by using the effective 
inverse slope parameter Tes instead of the temperature of the system together with the radial 
flow effect. The parameter 'jj^ also contains the density of spin and isospin states as well as 
the possible effect due to finite nucleon chemical potential. 

The inverse slope parameter in the nucleon pj, spectrum [Eq.(6)] is determined to be 
Teff = 295 MeV by using the experimental data of Ref. [32] for the minimum bias proton pj, 
spectrum. Taking the effective volume of the hadronic matter at freeze-out in minimum bias 
collisions to be \^ ~ 3, 700 fm'^, which is about 1/3 of the estimated hadronic matter volume 
in central Au+Au collisions at y/sj^ = 200 GeV [22], we obtain 7^ = 0.021 for the value 
of the parameter in Eq. (6).^ The comparison of our parameterized proton pj. spectrum, 
which is half of the nucleon pj, spectrum, with the experimental data is shown in panel (a) 
of Fig. 1. 

For the nucleon elliptic flow ^2^^ we assume that it is the same for both protons and 
neutrons and parameterize it in the form of the Fisher-Tippet (or Gumbel) distribution 
function [34], i.e., 

V2,n{Pt) =«7vexp{-exp[(Ajv-PT)//5iv]}- (7) 



The difference of the proton prp spectrum in minimum bias and in central collisions shows that the volume 
of the hadronic matter in the minimum bias colhsion is about 1/3 of that in the central collision when 
the data are parameterized using Eq. (6). 



6 




4 



1 2 3 4 1 2. 

Pj (GeV) (GeV) 



FIG. 2: (Color online) pj, spectrum (panel (a)) and V2 (panel (b)) of pions. Solid lines are 
parametrization of Eqs. (8) and (9). Experimental data are from Refs. [15, 32]. 



We use the experimental data of Ref. [16] for the proton V2, which were measured for 
20-60% centrahty, to fix the parameters. This gives = 0.258, = 0.683 GeV, and 
= 1.128 GeV. This parametrization is shown in panel (b) of Fig. 1 together with the 
experimental data of Ref. [16]. Our parametrization gives similar results in the measured prp 
region as the parametrization atanh(6pj. + c), which was introduced in Ref. [35] for quarks. 
But using the latter form would lead to a negative at small py. Although there is no 
experimental information on the nucleon V2 at very low pj, region, hydrodynamic model 
studies [36, 37] show that the nucleon t>2 produced in heavy ion collisions at RHIC has non- 
negative values [5]. Since the hydrodynamic calculations are expected to reasonably describe 
V2 at low prp region, we employ the form of Eq. (7) which gives non-negative V2,n even for 
very small prp. We may also take the proton V2 to have the form a' prp + j3'p^ for prp < p^ and 
ata.nh.{bprp + c) for prp > p^, which gives a similar result as Eq. (7) if we choose pj' = 0.75 
GeV, a = 0.293, b = 0.460 GeV-\ c = -0.175, a' = 0.0332 GeV~\ and p' = 0.0579 GeV'^. 



B. Pion spectrum and elliptic flow 

For the pion prp spectrum, the parametrization used for the nucleon prp spectrum [Eq.(6)] 
does not work well as it is affected by resonance decays [35, 38]. Instead, we use the 
parametrization, 

MPt) = 7. (1 + ^) [1 + V2,M cos(20)], (8) 

with = 1-29 GeV, /Jtt = —12.0, and 'Jt^ = 2.0 determined from the experimental data of 
the PHENIX Collaboration [32]. For the pion fg, it is parameterized as in Eq. (7), 

^2,^(Pt) = exp {- exp [(A^ - pp) (9) 

with = 0.184, = 0.461 GeV, and = 0.547 GeV to fit the PHENIX experimental 
data [15]. The comparison of our parametrization with experimental data is shown in Fig. 2. 
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VI. RESULTS 



With the above proton and pion pj, spectra and elhptic flows as inputs, we calculate in 
this Section the deuteron pj. spectrum and V2 from the dynamical model and compare the 
results with experimental data as well as those from the coalescence model. 



A. Deuteron spectrum and elliptic flow from two-body reactions 

In terms of rapidity and transverse momentum, Eq. (3) for the deuteron production rate 
can be written as 

d^TZ 1 f ^ 

d?p^^dy^ " 8(2^ J \{^y^^PT,id^iPT,ifN{.PTM + UP^+P^-Pd)\ 

X \M{NN ^ d'K)\^5{{p,+p2-p,f - Ml). (10) 

The deuteron pj, spectrum is then obtained by multiplying the above rate by the volume 
of the hadronic matter at freeze out V ~ 3,700 fm^ and the time interval At ~ 4 fm/c. 
The latter value is chosen to reproduce reasonably measured deuteron spectrum. It takes 
into account the fact that in a more microscopic approach that follows the production and 
dissociation of deuterons from the reaction A^A^ — > as well as the deuteron elastic scat- 
tering, deuterons are expected to be produced over a finite time interval from the expanding 
hadronic matter formed in heavy ion collisions. To extract the elliptic fiow V2 of produced 
deuterons, we further express the result as 

fd{PT,d^ yd)[l + '2v2,d{PT,d) cos(2$<i)], (11) 



PT,d^PT,d^^ddya 

where is the azimuthal angle of produced deuteron and 0^ 2 those of the two nucleons 
relative to the reaction plane. Because of the 5-function which comes from the energy- 
momentum conservation and the pion on-shell condition, one of the nucleon azimuthal angles 
is constrained by kinematics. This is different from the assumption in the naive coalescence 
model that the direction of the nucleon momentum is constrained to (pi = (f)2 = ^d- 

The resulting deuteron pj, spectrum using the phenomenological scattering amplitude 
given in Appendix A for the NN — > dir reaction is shown by the solid line in panel (a) 
of Fig. 3. Compared with the measured spectrum of deuterons shown by filled circles, 
ours is softer. This result indicates that the radial fiow effect on produced deuterons is not 
fully taken into account in the dynamical model. Since the cross sections for ird elastic 
and inelastic scatterings are comparable at pion kinetic energy of around 100 MeV [39], 
deuterons produced from the dynamic model are expected to suffer further scatterings and 
thus to approach thermal equilibrium with pions. To properly treat the radial flow effect 
on produced deuteron spectrum requires the inclusion of the processes considered in present 
study in a non-equilibrium transport model, such as the ART model [40]. As we are mainly 
interested in studying the effect of energy conservation on the deuteron elliptic flow, we will 
leave such more realistic approach to a future study. Surprisingly, the deuteron px spectrum 
obtained from Eq. (1) of the coalescence model agrees better with the data as shown by the 
dashed line. 
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FIG. 3: (Color online) (a) pj< spectrum of deuterons. Experimental data are from Ref. [33] for the 
minimum-bias collisions (filled circle), (b) Scaled deuteron elliptic flow V2 as a function of prp/2 at 
= 0. The dotted line is the nucleon elliptic flow V2 j\f of Eq. (7). Experimental data are from 
the PHENIX Collaboration [16] (filled circles) and the STAR Collaboration [17] (filled squares). 
In both figures, the dashed and solid lines are results from the coalescence model of Eq. (1) and 
the dynamical model in the present work, respectively. 



Results from the dynamical model [Eq. (1)] for the scaled deuteron V2, i.e., ^2^/2, are 
shown by the solid line in panel (b) of Fig. 3 as a function of Compared to the exper- 

imental data shown by filled circles from the PHENIX Collaboration and by filled squares 
from the STAR Collaboration, our model can describe both reasonably well, particularly 
those from the PHENIX Collaboration [16]. The calculated deuteron is, however, always 
positive, which is not consistent with the negative values of deuteron V2 at small prp seen 
in the preliminary STAR data [17]. This may again be related to our neglect of final-state 
interactions of produced deuterons. If the nucleon number scaling of deuteron f 2 is exact 
as in the naive coalescence model, the scaled deuteron elliptic flow f2d(PT/2)/2 should be 
identical to the nucleon elliptic flow V2n{Pt) \^'^- i^)] shown by the dotted line. Although 
some deviations are found between the two curves at large pj. region, our result supports 
the idea of the nucleon number scaling of V2 a- Furthermore, our result is very similar to 
that of the more general coalescence model for deuteron production (dashed lines). In the 
considered energy region, the two approaches give nearly the same f 2 ^. This may indicate 
that the momentum conservation in the coalescence model plays a more important role than 
the energy conservation that is not respected by the coalescence model. Therefore, it would 
be important to clarify the differences between our dynamical approach and the coalescence 
model. 

As we have discussed before, even in the more general coalescence model [Eq.(l)], the 
momenta of the constituent particles are still nearly aligned to the momentum of the pro- 
duced particle. In the case of deuteron, this is because the deuteron wave function prefers 
nearly equal momenta for its constituent nucleons, i.e. = in Eq. (2). Combined with the 
momentum conservation, this leads to ~ ~ P25 but keeping = +P2- Therefore, 
although employing the physical deuteron wave function can relax the momentum constraint 
p^/2 = Pi = P21 the momenta of the constituent nucleons remain restricted and take only 
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FIG. 4: (Color online) The allowed regions of and P2 in A^A^ dir reaction. The empty space 
represents the unphysical phase space. Left panels show the phase space in the azimuthal angles 
and <j)2 of the two nucleons and right panels show the phase space in \prp^\ and |P2"2l- Upper 
panels are for = 0.1 GeV and lower panels for Prp^ = 4.0 GeV. The deuteron azimuthal angle 
is fixed at = 0. 

values around the region that appears in the naive coalescence model, which in turn brings 
the prp dependence of deuteron t>2 to scale approximately according to its constituent nucleon 
number. 

In the dynamical model, the role of the deuteron wave function appears through the 
squared transition amplitude for A^A^ — dir, where the energy is conserved. As a result, 
the allowed momenta of the nucleons forming the outgoing deuteron are different from those 
in the coalescence model. This can be seen in Fig. 4, which shows the allowed momenta 
Pj^i and pj^2 of nucleons to form a deuteron of pT,d = 0.1 GeV (upper panels) and 4.0 GeV 
(lower panels) with $d = 0. The naive coalescence model allows only one point in each 
figure, namely 0i = 02 = in the left panels and prpi = Pt2 = VtaI'^ ^^e right panels. 
The more general coalescence model would allow a larger region, but the major contribution 
comes from the region of {Pti^Ptq) ~ KVt J'^iVt J'A (0i,02) ~ (0,0) (modulo 27r). 
The dynamical model based on the A^A^ — > (ivr reaction covers, on the other hand, the 
phase space shown by the shaded areas in the figure, which may or may not include those 
involved in the coalescence model. For the phase space in transverse momenta, because of 
the exponential form of the nucleon prp distribution function, the main contribution in the 
dynamical approach still comes from that allowed in the more general coalescence model, 
i.e. Pxi ^ Pt2- For the phase space in the azimuthal angles, one can find that the region 
allowed in the dynamical model for pT,d = 4.0 GeV (the lower-left panel of Fig. 4) covers 
that of the coalescence model, which gives the major contribution to deuteron production. 
This explains the resemblance of our result with the coalescence model in the intermediate 
and large pj, regions. 

The similarity between the two models does not hold, however, in the low pj. region. As 
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shown in upper panels of Fig. 4 for = 0.1 GeV, the kinematically allowed phase space 
in the azimuthal angles is very different from the allowed region in the coalescence model. 
In fact, the physical phase space excludes the region of (0i,02) ~ (0)0) (modulo 27i) as a 
result of energy- momentum conservation and the phase space of the coalescence model is not 
physically allowed. Therefore, the energy conservation, which is neglected in the coalescence 
model, plays an important role in the low pj, region as one would have expected. In the 
dynamical model based on the NN — > di: reaction in heavy-ion collisions, there is thus no 
a priori reason to expect to have the constituent number scaling of the deuteron V2 at low 
prp region. 

The scaling behavior of the deuteron V2 shown in the low region in panel (b) of Fig. 3 
is a coincidence coming from the parameterized nucleon V2 used in the calculation. In 
fact, if we allow negative at very low pj. region, the coalescence model would predict 
a negative V2 ^ as expected from the scaling behavior. However, in the present dynamic 
model calculation, although inclusion of deuteron final-state interactions is still required, 
this expectation is not valid any more in the low pj, region, and we always have small but 
positive V2 d in this region even if we assume negative nucleon V2 for low pj,. This is because 
of the fact that the nucleon azimuthal distribution function is weighted by the nucleon 
and its distribution function as can be seen in Eq. (10). Therefore, for p^,^ = 0.1 GeV, the 
main contribution to deuteron production comes from nucleons of pj, ~ 0.5 GeV, while it 
must be ~ 0.05 GeV in the coalescence model due to the constraint on the magnitude of 
the momenta. Furthermore, the physical azimuthal angles in the dynamical model do not 
include the region allowed in the coalescence model. Consequently, these cause the main 
difference between the coalescence model and the dynamical model in the present work. 

B. Deuteron spectrum and elliptic flow from three-body reactions 

For deuteron production from the three-body NNN — > dN and NNtt — > dTi reactions, 
the deuteron transverse momentum spectrum and differential elliptic flow can be similarly 
calculated as in the case of the two-body NN — > dn reaction. Multiplying the calculated 
rates by the same hadronic matter volume K ~ 3, 700 fm^ at freeze out and the time interval 
Ar ^ 4 fm/c, the results are shown in Fig. 5 by the dashed and dotted lines, respectively. 
For comparison, results based on the two-body NN — > d7r reaction, which have been shown 
in Fig. 3, are also given by solid lines. 

Although the physically allowed nucleon phase spaces for deuteron production in NNN 
and NNtt reactions are different from that in the A^A^ dTi reaction, the predicted deuteron 
f 2 from the three reactions do not show drastic difference from each other as shown in panel 
(b) of Fig. 5, and all lie close to the experimental data except at very low momentum where 
the deuteron V2 is negative in the preliminary data from the STAR Collaboration [17] . 

In contrast to the deuteron V2, the predicted deuteron transverse momentum spectrum 
depends strongly on the reaction mechanism. As shown in panel (a) of Fig. 5, contributions 
from three-body reactions are suppressed by several orders of magnitude compared to that 
from the two-body reaction. This is mainly due to the low nucleon and pion densities at 
freeze out. Another possible, but subsidiary, reason for this behavior is the magnitude of 
the reaction amplitudes used in our study. The very simple model we have used for the 
NNN — > dN and NNtt — > dn reactions may have underestimated their contributions to 
deuteron production. On the contrary to the phase space in NN — > dn reaction (Fig. 4), 
the phase space in the three-to-two-body reactions has no restrictions on the angles and 



11 




FIG. 5: (Color online) Calculated deuteron pj, spectrum (panel (a)) and scaled V2 (panel (b)) at 
mid-rapidity from A^A^ — > dn (solid lines), NNN — > dN (dashed lines), and NNtt — > dn (dotted 
lines) reactions. Experimental data are from Refs. [16, 17]. 



magnitudes of the momenta of three initial particles. 



VII. SUMMARY AND DISCUSSIONS 



In this paper, we have studied the deuteron elliptic flow in heavy-ion collisions at RHIC 
using a dynamical model for deuteron production that conserves both momentum and en- 
ergy. Since elliptic flows in relativistic heavy ion collisions are mainly generated during early 
partonic stage, deuteron elliptic flow is expected to not change much during the hadronic 
stage, we have considered the deuteron production rate from the NN —>■ dn reaction at 
the thermal freeze out. The deuteron yield is then obtained from multiplying the produc- 
tion rate by the volume of the hadronic matter at freeze out and a time interval of about 
4 fm/c. Because of the neglect of deuteron final-state interactions, which would enhance 
the incomplete treatment of radial flow effect, the resulting deuteron pj, spectrum is some- 
what softer than that measured in experiments. The obtained result for the deuteron V2 
agrees, however, reasonably with the measured one, except at very low where the latter 
is negative. Our results further show that this model can reproduce the predictions from the 
coalescence model. The origin of this result could be understood from the physically allowed 
phase space of nucleon momenta in deuteron production. At large pj., the phase space in 
our dynamical model covers that of the coalescence model, which eventually dominate the 
production process because of the nucleon distribution functions. Therefore, our dynamical 
model supports the description of the coalescence model. However, at very low pj,, the 
nucleon phase space in the coalescence model turns out to be physically forbidden because 
of energy-momentum conservation. Thus, deuteron production via NN — ^ dn generally 
does not reproduce the prediction of the coalescence model at very low prp. Although the 
prediction from the dynamical model follows the scaling behavior, this property is strongly 
dependent of the nucleon f 2 at the low pj, region. 

We have also considered three-body NNN dN and NNtt —>■ rfvr reactions for deuteron 
production in heavy-ion collisions. Because of the different nucleon phase space involved 
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N(p2) N(q) d(d) 



FIG. 6: Diagram for the NN dir reaction. The blob denotes the ttN ttN elastic scattering. 

in these reactions from that in the two-body NN — > dn reaction, these reactions lead to 
different predictions on the deuteron The deuteron production rates from these reactions 
are, however, suppressed compared to that from the two-body A^A^ — > dn reaction by several 
orders of magnitude. This implies that deuteron production in heavy-ion collisions is mainly 
due to the NN dn reaction, which has the largest cross section among the NN — > dX 
reactions [23]. We have also found that all models considered in the present study give 
positive deuteron at very low p^. This is in contrast to preliminary data from the 
STAR Collaboration that deuterons produced in heavy ion collisions at RHIC may have 
negative V2 at very low pj, region. To obtain negative deuteron V2 required negative nucleon 
V2 in the coalescence model. Therefore, precise measurements of both the nucleon and 
deuteron f 2 in this transverse momentum region are highly desirable for understanding the 
deuteron transverse momentum and elliptic flow in heavy-ion collisions. For our dynamical 
model approach, there are several issues which deserve further studies. As mentioned in 
Section VI, final-state interactions of deuterons, particularly the nd elastic scattering, need 
to be taken into account via a microscopic transport model or a schematic kinetic model. 
Also, the deuteron production amplitude used in the present study assumes a point-like 
structure of the deuteron. For a more reahstic description, the spread of the deuteron wave 
function should be taken into account in the production amplitude to render a more direct 
comparison with the coalescence model. 
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APPENDIX A: NN dn REACTION 

The A^A^ dTi reaction can be described by the one-pion exchange diagram shown in 
Fig. 6. In this model, one of the incoming nucleons emits a virtual pion which causes the 
ttN — > ttN reaction shown by the blob in Fig. 6 before forming the outgoing deuteron. 
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There are other diagrams which include the nucleon exchange diagrams but it is well-known 
that their contributions are negligible [41]. In this study, we thus adopt a simple model for 
the NN — > dn reaction of Fig. 6 following the method of Refs. [42, 43]. Here we focus on 
the reaction of pp dir'^ and use the isospin relations a{pp — > dir'^) — a{nn — > dir') — 
2a{pn — > dn^) to obtain the isospin averaged cross section a{NN — >• dn) = a{pp dn^). 

Using the momenta indicated in Fig. 6, the amplitude for the pp dTr"*" reaction can be 
written as 



1 -g ■ 7 + 



dNN 



N 



{d-q)-j + M, 



X7lf^peif(^ + S*-7)«W. (Al) 

where Mn and Mj^ arc, respectively, the nucleon and pion masses, and FTrA^Ar and TdNN 
denote, respectively, the pion-nucleon and deuteron-nucleon vertices. The functions A and 
B depend on the Mandelstam variables and are related to the invariant nN elastic scattering 
amplitude by 

MinN TiN) = u{pf) {A + Bk- 7) m(pJ, (A2) 

where p^ and pj are the initial and final nucleon momenta, respectively. The squared matrix 
element can be written in terms of the ttN elastic scattering cross section as 

\M{7rN nN)\' = (87r)Xiv^(5.iv, «.^)- (A3) 

where the Mandelstam variables for the vrA^ system are s^jy = (^j+Pi)^ and u^j^ = {k^—pjY, 
with fcj being the momentum of the initial pion. In our qualitative study, the elastic nN 
scattering can be well approximated by assuming that it is dominated by the A resonance 
with the interaction Lagrangian 

>C.ArA = ^A^d.nN + H.c, (A4) 

where fnNA = 2.23 and the isospin factors can be found, e.g., in Ref. [44]. 
Following Ref. [41], the spin averaged squared amplitude is then given by 

\M\' = 

4 ^ ■ • 

spins 

2") 

\M{TrN ^TTN)f. (A5) 

In the above, g (~ 4%) is the amount of D-wsve state in the deuteron wave function, 
is the deuteron mass, and k = \/Mj^e with the deuteron binding energy e (r^ 2.23 MeV). 
The pion-nucleon coupling constant is denoted by Qt^nn with a value glj^j^/ii: — 14.0, and 
GTrNN{(f) is the form factor 

_ , - ^ , ^ 
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FIG. 7: (Color online) Total cross section for the pp — > dn^ reaction. The experimental data are 
from Refs. [45-48]. 



with A = 0.73 GeV. The deuteron- nucleoli coupling constant g^i^i^ is given by 



9dNN 



[MN{l-nK){l + Q^) 



1/2 



(A7) 



where r^ 1.75 fm) is the effective range for the spin-triplet nucleon-nucleon elastic scat- 
tering. To reproduce the experimental total cross section for pp dn^ within this simple 
model, we use Ma = 1200 MeV. The calculated cross section is shown in Fig. 7 and is seen 
to reproduce very well the measured one. 



APPENDIX B: NNN dN AND NNtt dn REACTIONS 

The amplitudes for producing deuterons from three-body reactions NNN dN and 
NNtt dir are approximated by making use of the amplitude of NN — » dn obtained 
above. For the NNN — >■ dN reaction, we consider the one-pion exchange process, namely, 
the collision of two initial nucleons gives the deuteron and an off-shell pion, which is then 
absorbed by the spectator nucleon. In this approximation, the reaction amplitude becomes 

\M{NNN dN)\' = {p,-Pf - Ml) \M[NN d^)\\ (Bl) 

where is the intermediate pion momentum, is the momentum of the initial spectator 
nucleon, and p ^ is the nucleon momentum in the final state. 

Replacing the off-shell pion by an off-shell rho-meson, which is then absorbed by a specta- 
tor pion leads to an approximate description of the NNn dir reaction. The amplitude for 
the subprocess NN dp can be obtained by the same method of previous section for obtain- 
ing the A^A^ —>■ dn amplitude. The result is thus similar to Eq. (A5) with |A^(7rA^ — > 7rA^)p 
replaced by {AilnN —>■ pN)\'^/2. A complete analysis of the vrA^ — >■ pN reaction requires 
many possible terms arising from meson exchanges and nucleon resonance contributions. 
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FIG. 8: (Color online) Total cross sections for vr+p p^p (panel (a) and pp dp'^ (panel (b)). 
Experimental data are from Ref. [49]. 



For a simple estimate on the cross section for this reaction, we describe it by one-pion ex- 
change only and fit the experimental data with suitable form factors. The relevant effective 
Lagrangians are 



r 

^ttNN 



9ttNN 



2Af, 



N 



(B2) 



where Qp^^ = 6.02 and ^^^^/(47r) = 14.0. The form factor in the form of A7[A^ + (t- M^) 2] 
is used with A = 0.65 GeV. The NNn —>■ dn amplitude is then given by 



\M{NNn ci7r)p 



ill - M^f 



(B3) 



Our results on the total cross sections for vr+p p^p and pp — > dp^ are shown in Fig. 8. 
The former is seen to agree with the measured cross section. 
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